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ABSTRACT. A far-wing line shape theory based on the binary 
collision and quasistatic approximations that is applicable 
for both the low- and high-frequency wings of allowed 
vibration-rotational lines has been developed. This theory 
has been applied in order to calculate the frequency and 
temperature dependence of the continuous absorption 
coefficient for frequencies up to 10 000 cm -1 for pure H 2 0 
and for H 2 0-N 2 mixtures. The calculations are made assuming 
an interaction potential consisting of an isotropic Lennard- 
Jones part and the leading long-range anisotropic part, and 
utilizing the measured line strengths and transition 
frequencies. The results compare well with existing data, 
both in magnitude and in temperature dependence. This leads 
us to the conclusion that although dimer and collision- 
induced absorptions are present, the primary mechanism 
responsible for the observed water continuum is the far-wing 
absorption of allowed lines. Recent progress on near-wing 
corrections to the theory and validations with recent 
laboratory measurements are discussed briefly. 


1* INTRODUCTION 


The absorption of radiation by water vapor in the window 
regions of the Earth's atmosphere is a long-standing problem 
[1] of great practical as well as theoretical interest [2]. 
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In the past few years there have been a number of general 
review [3-6] concentrating primarily on the experimental 
characterization of this absorption, which for historical 
reasons is called the "water continuum", and the interested 
reader is referred to these reviews for original literature 
citations. As a result of this extensive work, there is 
nearly unanimous agreement on the density dependence 
(quadratic) , general agreement on the temperature dependence 
(strong, negative, at least in the 8-13 fm window for 
atmospheric temperatures) , but considerable disagreement as 
to the magnitude and physical mechanism responsible for the 
absorption. Water dimers [7] and polymers [8], collision- 
induced absorption [ 9 ] , and the superposition of the far 
wings of collisionally broadened allowed dipole lines [1] 
have all been suggested as possible sources of the water 
continuum. While the first two mechanisms are undoubtedly 
present and contribute to the absorption in the atmosphere, 
there is increasing evidence that they are not the primary 
mechanism. 

Within the past decade, substantial progress in 
calculating the absorption in the far wings of allowed 
transitions has been made. Starting from the basic 
formalism of Fano [10] and Davies et al. [11], Rosenkranz 
[12,13] has developed a quasistatic statistical theory that 
is applicable in the high-frequency wing of the pure 
rotational band (the 8-13 jon window) . This theory 
successfully predicted both the magnitude and the strong, 
negative temperature dependence. A line-by-line 
generalization of Rosenkranz ' theory was proposed by Boulet 
et al. [14], but because of the extensive computations 
involved, it has been implemented only for atom-molecule 
collisions. In a series of recent papers [15-20], Ma and 
Tipping have eliminated many of the restrictions in the 
quasistatic theory so that it now constitutes a unified 
theoretical framework within which one can calculate the 
magnitude and temperature dependence of the self- and 
foreign- broadened far-wing absorption. 

In the present paper, we briefly review the 
mathematical formalism and the nature of the approximations 
made. We then discuss the theoretical line shapes and the 
resulting absorption coefficient a(o) for the case of pure 
water vapor, and present comparisons with selected 
experimental data. Finally, we indicate the principal 
difficulties of the present theory, and discuss possible 
improvements that are feasible within the current framework. 


2. THE GENERAL THEORETICAL FRAMEWORK 


For a low-density gas in thermal equilibrium at a 



temperature T, binary collisions dominate and one can focus 
on one absorber molecule and treat the other molecules as a 
bath. In this case, the absorption coefficient per unit 
volume a(u) at frequency cj is given by the well-known 
formula 


a(u) = (4n 2 /3hc) n a o tanh(fuj/2kT) 

x [F(u) + F (-«) ] , (1) 

where the spectral density, F(u), is given by 

F(u) = tt" 1 Re Tr | e lut <J a (0) *p a (t)> dt . (2) 

J 0 

•4 

In this expression, p a (t) is the dipole operator of the 
absorber molecule in the Heisenberg representation and the 
angular brackets denote the ensemble average over all (one 
absorber plus all bath molecules) variables. Introducing 
the total Liouville operator L 

L * Lq + = + L^) + , (3) 

where L 0 is the unperturbed part describing the absorber a 
and bath molecules b, and L, corresponds to the interaction 
between them, the spectral density can be written in the 
form 


F(w) - i Re Trj% lwt {J a (0)-e _lLt p p a (0)}dt , (4) 

where p is the density operator. As shown in [18], e _iLt 
can be approximated by 

e -iLt = e -iL 0 t e ”i L it 


where 


x {1 



L' 1 (t 1 ) + *••} , 


(5) 


L'(t) = e iL l t e iL 0 t L x e" iL 0 t e _iL l t - L x . (6) 

Substituting Eq. (5) in Eq. (4) , the spectral density can be 
represented by the expansion 

F(u) = F q (<j) + F^u) + ... 

where F^o) is the first-order correction term resulting 


( 7 ) 



from the non-commutation of L 0 and L^. Within the 
quasistatic framework, F 0 («) and F 1 (<j) can be expressed in 
termg of symmetrized line-shape functions [15-20] Xi-j ” ( u ) 
and respectively, where the transitions are 

labeled by subscripts i and j; in principle, these can be 
calculated from the interaction potential V(?) if we assume 
it to be of the form 

v (r) = V iso (r) + V aniso (r) , (8) 

and only consider the leading long-range part of the 
anisotropic potential (dipole-dipole for H 2 0-H 2 0 or dipole- 
quadrupole for H 2 0-N 2 ) . This approximation simplifies the 
theoretical expressions for the line shape functions but 
limits the accuracy of the results. Recently, we have 
generalized the theory to incorporate multi-component 
potentials, but have not yet carried out calculations for 
H 2 0 [20]. To reduce the computation time required, one 
usually introduces band-average line shape functions 
X' 0 ' (cj) and x 1 («) > that is, the same line shape is assumed 
for every line in a given vibrational band. Because the 
derivations as well as the theoretical expressions for the 
individual and band-averaged symmetrized line shape 
functions are rather complicated, they are not given here 
but can be found in [15-20]. 

Using the results for the spectral densities expressed 
in terms of the symmetrized line shape functions, one can 
write the absorption coefficient a(u) in the form 


where 


u sinh(ha/2kT) 


a(u) = n a 2 S ^ 4 . : — - — 

cj^j>0 D u ij sinh(h« i j/2kT) 


x n ~ 1 {■ 


(a - Uij) 


2 Xf(« “ «ij) 


(w + Wij) 


- X-(« + «ij) ) # 


X+(") = X+ (0) (") + X+ (1) («) 


(9) 


( 10 ) 


and S^4 are the usual line strengths. The subscripts ± 
indicate the positive- or negative-resonance line shape 
functions [15-17]; these have a different functional 
dependence on the frequency and reflect the fact that 
neither the individual nor the band average line shape 
functions are symmetric functions. If the differences 



between the positive- and negative-resonance line shape 
functions are neglected, one gets the band average line 
shape function introduced by Rosenkranz [12,13], which is 
written without a subscript. 

The results for x'°* («) , X^ 1 ' («) / and x(«) are shown in 
Fig. 1 for pure water calculated at T = 296 K. If the lines 
were Lorentzian, these curves would be horizontal lines 
whose intercepts would be the corresponding half-widths. 
However, as can be seen from the figure, the shapes are 
super-Lorentzian out to approximately 200 cm" 1 from the line 
center, after which they fall off very rapidly with 
frequency. 



Figure 1. The line shape functions (in units of cm” 1 atm" 1 ) 
as a function of frequency u (in units of cm" 1 ) calculated 
at T * 296 K. The dotted curve is the band-average line 
shape function without the near-wing correction, the dashed 
curve is the near-wing correction, and the solid 
curve is the total. 


Furthermore, the correction term is significant between 10 
and 80 cm” 1 , after which it is negligible. For 
displacements less than 10 cm” 1 it is comparable to or even 
greater than the leading term, but in this region the 
quasistatic approximation used to calculate the line shape 
functions is not valid [19]. 

By excluding any contribution from frequency 
displacements less than 25 cm” 1 , one can calculate the 
absorption coefficient using Eq. (9) ; the results are shown 
in Fig. 2 for T = 296 K. One can clearly see that the 
magnitude of the far-wing contribution to the absorption 
varies over 6 decades in the frequency region 0-10 000 cm” 1 . 



FREQUENCY 


Figure 2. The absorption coefficient a(u) (in units of cm 2 
molecule” 1 atm” 1 ) as a function of frequency u (in units of 
cm” 1 ) calculated at T = 296 K; the solid line is the result 
using two symmeterized line shape functions, while the 
dotted line is the result using the average line shape 
function. 



3. COMPARISONS BETWEEN THEORY AND EXPERIMENT 


In order to illustrate the accuracy of the 
calculations, we present some comparisons between theory and 
experiment; more comparisons are given in [15-19]. In Fig. 

3 we show results for the 300 cm -1 < u < 1100 cm -1 region at 
T = 296 K. Similar results obtain for other temperatures in 
this region and for other spectral regions, although the 
agreement is not as good as that shown. However, in these 
other spectral regions, the continuum contribution is much 
weaker and the corresponding experimental uncertainties are 
accordingly larger. 
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Figure 3. The absorption coefficient a(u) (in units of cm 2 
molecule" 1 atm" 1 ) as a function of frequency for 300 cm" 1 
< cj < 1100 cm" 1 at T = 296 K. The experimental values of 
Burch et al. [21] are denoted by +, while o and ▲ correspond 
to the theoretical results obtained without and with the 
near-wing corrections , respectively. 



Furthermore, to obtain the experimental continuum, one has 
to subtract out the contributions from nearby lines, and 
different authors use different methods (i.e., different 
cutoff displacement and/or different line shape functions to 
represent the central parts of the lines) . Nevertheless, as 
can be inferred from these comparisons, there is generally 
good agreement in the magnitude of the absorption. 

In Fig. 4 we compare the measured and observed 
temperature dependence. Although the magnitude of the 
theoretical continuum is less than that observed, the shape 
is in good agreement. 
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We note that for atmospheric temperatures, there is a 
strong, negative temperature dependence, but this dependence 
weakens as T increases so that for T > 400 K, the dependence 
becomes weakly positive. This latter behavior is in sharp 
disagreement with the exponential decrease one would expect 
from dimer absorption. 

In the comparisons shown above and those not shown but 
published previously in [15-19], there are no adjustable 
parameters except the Lennard- Jones parameters describing 
the isotropic part of the interaction potential. These can, 
in principle, be determined from other types of 
measurements, but in our previous calculations, we have 
adopted the same potential as reported by Rosenkranz [12,13] 

v iso ( r ) * c ff ‘ 6 C(*/r) 48 - (*/r) 6 ] (U) 

where C/k * 8 x 10 5 X 6 K, and a m 3.1§ A. Recently we have 
carried out calculations in the 700 cm" 1 to 1300 cm" 1 
frequency region for comparison with accurate experimental 
data of Kulp [25]. In Fig. 5 we plot the experimental 
transmission spectrum along with a line-by-line synthesis 
including the theoretical continuum absorption that is shown 
separately as well. This was calculated using a 6-12 
Lennard-Jones isotropic potential 

v iso ( r ) * 46 [(*/*) 12 - (*/r) 6 ] (12) 

with e/k = 600 K and a - 3.12 A; this potential is about 8 % 
deeper than that in Eq. (11) and significantly wider. 

From the difference between the experimental data and the 
model plotted in the lower part of the figure, one can see 
that there is no systematic difference and the overall 
agreement is very good. 


4. DISCUSSION AND CONCLUSIONS 


The results discussed in the previous section for the 
line shape functions and continuum absorption coefficient 
(together with more extensive results given in [15-19]) 
clearly show that the far-wing absorption resulting from an 
interacting pair of water molecules cam account for both the 
magnitude and the temperature dependence of the observed 
continuum. This absorption has been calculated by making 
the following approximations. 1) Only binary collision are 
considered; for low density gases, this approximation is 
certainly valid in the wings of the lines. 2) The 
quasistatic approximation is used; by introducing the 


Liouville formalism, the Fourier transform in the spectral 
density is carried out formally and time does not appear 
explicitly. The trace over the translational degrees of 
freedom is replaced by an integration over the separation r 
between the centers of mass of the colliding pair using the 
pair potential as a weighting function; the integrand is 
the matrix elements of the resolvent operator pertaining to 
the anisotropic interaction which must be diagonalized. 
Because of this, only one term with a specific spherical 
symmetry was treated; if more than one term is important, 
other nonessential approximations can be introduced in order 
to make the problem tractable [20]. 


MICRONS 



Figure 5. The water vapor transmission spectrum in the 700 
cm” 1 - 1300 cm” 1 spectral region measured at a pressure of 
8.1 torr at T = 297 K using a path length of 252 m [25]. 

The dotted curve represents the measurements while the solid 
curve is a synthetic spectrum calculated by a line-by-line 
method using 92 HITRAN molecular data [26] and the 
theoretically calculated continuum which is shown 
separately. The difference in absorbance between the model 
and the data is plotted in the lower part of the figure. 




3) Other approximations such as the band average 
approximation have been introduced for computational 
simplicity? these could be avoided if need be, but in light 
of the other approximations and the extra computation time 
required, this is not warranted at present. 4) As indicated 
above, we have only considered simple Lennard-Jones models 
for the isotropic potential. Other more accurate models or 
even numerical potentials could be easily incorporated 
within the present framework. 5) When comparing calculated 
spectra with experiment (cf. Fig. 5) the effects of the 
local lines must of course be included. To do this, one has 
to assume some profile for the absorption near the line 
center and extrapolate to the near wings (say to 10 or even 
25 cm" 1 ) . In this intermediate region, neither the impact 
nor the quasistatic theory is applicable and the true line 
shape is not known. By including the near-wing effects 
represented by the spectral density term F^u), this 
extrapolation is reduced although the basic problem remains. 
For the window region around 1000 cm" 1 this problem is not 
crucial, but closer to a band center, it can be important 
[19]. 

In conclusion, we note that much progress on the 
theoretical calculation of the far wings of allowed lines 
has been made. The far-wing absorption is the primary 
mechanism responsible for the water continuum and this is 
described reasonably well within the framework of the 
present theory. However, improvements to the theory in 
order to calculate the absorption due to the near wings of 
spectral lines, the use of multi-component anisotropic 
potentials [20] and more accurate isotropic potential 
models, and the possible inclusion of other mechanisms (for 
instance, interferences effects between allowed and induced 
dipole moments) are needed. Finally, we would like to 
stress the need for more accurate experimental data 
(laboratory and atmospheric) that can be used to assess the 
accuracy of the theoretical results. 
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